Gouy phase for full-aperture spherical and cylindrical waves 
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We investigate the Gouy phase shift for full-aperture waves converging to a focal point from all directions in 
two and three dimensions. We find a simple interpretation for the Gouy phase in this situation and show that 
it has a dramatic effect on reshaping sharply localized pulses. © 2012 Optical Society of America 
OCIS codes: 070.7345, 050.5080 



Gouy phase shift [1,2] has been known for more than 
a century and it still attracts a lot of interest in the opti- 
cal community. Its origin has been explained in different 
ways and contexts, for an overview see Ref. [3] and the 
references therein. 

Gouy phase is a phase shift of a converging wave ob- 
tained for instance by focusing light with a lens when 
passing through the focus. It turns out that the phase 
change of the wave upon passing through the focal point 
is smaller by it compared to the situation if a plane wave 
were propagating instead. In other words, the local wave- 
length near the focus is slightly larger than Ao = 2itc/uj, 
the wavelength of a plane wave of the same frequency 
uj. One interpretation of this fact is that near the fo- 
cus the wavevector has non-negligible transversal com- 
ponents and therefore its longitudinal component must 
be somewhat reduced, which increases the local wave- 
length [4]. 

Gouy phase is usually discussed for light that propa- 
gates more or less in one direction, as when focusing light 
with a lens or creating a Gaussian beam. On the other 
hand, recently other geometries of light propagation have 
been studied in relation to perfect imaging devices such 
as Maxwell's fish eye [5,6] or to Air microscopy [7]. There, 
light arrives at the focal point from all spatial directions 
and/or leaves that point for all directions. In this pa- 
per we show that Gouy phase exists in this situation 
as well and has a nice and straightforward interpreta- 
tion. We also discuss the influence of Gouy phase shift 
on light pulses passing through the focal point and show 
that it has a dramatic effect on the pulse shape in the 
two-dimensional case. 

Our starting point is the wave equation 

c 2 Avb - fa = (1) 

that can be used for describing different types of waves, 
e.g. scalar or electromagnetic waves, in a homogeneous 
non-dispersive medium. Here c is the speed of the waves. 
Now consider a monochromatic spherical wave of fre- 
quency uj converging to the focal point at the origin of 
coordinates. This wave can be expressed as 



with k — ui/c. However, Eq. (2) is not a solution of 



Eq. (1) at the very origin because the left-hand side di- 
verges there. Therefore, if there is no absorber (drain) 
for the radiation at r = 0, Eq. (2) does not represent 
the full solution of the wave equation. This has a simple 
reason: the converging wave cannot just disappear at the 
focal point [6] but it changes there into a diverging wave 
Vwt — &cxp[ifcr — \u)t]/r which must be superimposed 
with (2) to get the full solution. The only way to satisfy 
the wave equation at the origin is to set b — —a, i.e., add 
the diverging wave with a phase shift 7r, which gives the 
total wave 

V- 1 — ipin + Vwt = — 2ifcasinc (kr) exp(— iut) (3) 

The phase shift of it between the two waves at the origin 
can be interpreted as the Gouy phase. 




Fig. 1. Comparison of the function sine /ex (blue thick 
curve) and — sin kx (red thin curve) describing a super- 
position of converging and diverging spherical waves and 
a plane wave, respectively, along an axis passing through 
the focal point. The phases match in the region x < 
while they differ by ir for x > 0, which demonstrates the 
Gouy phase shift. 

To illustrate this in another way, let us have a look 
more closely at the wave in Eq. (3). Fig. 1 shows two 
functions: one is sine kx, which represents the spher- 
ical wave in Eq. (3) (apart from the global factor 
— 2ifca exp (— iwt)) along an x-axis passing through the 
focal point. The second function is (— sin /ex) and rep- 
resents a plane wave of the same frequency propagating 
along the x-axis whose phase was chosen to coincide with 
the phase of the sine wave in the region x < 0. Obviously, 
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in the region x > the phases of the two waves differ: 
the sine wave is delayed by 7r with respect to the plane 
wave. We see that a similar thing happens here as near 
the focus of a lens. This demonstrates the Gouy phase for 
spherical waves in an elegant way, by simply comparing 
the graphs of the functions sine kx and sin kx. We can 
also see that the phase shift occurs on the length scale of 
the order of the wavelength, which is consistent with the 
case of the lens where the Gouy phase shift occurs on the 
scale of A/ 2 /a 2 with / denoting the focal distance and a 
the aperture [8]; in our case / and a can be considered 
to be of the same order. 




Fig. 2. Similar as Fig. 1, but for the 2D situation. The 
blue thick curve shows the cylindrical wave Jo(fcx) and 
the red thin curve shows the plane wave sin(fcx + 3tt/A). 
The phases match for x ;$ — A = — 2-k jk but differ by 7r/2 
for x j§ A - a clear demonstration of the Gouy phase shift 
in two dimensions. 

A similar consideration can be made for the two- 
dimensional rotationally-symmetric waves (or cylindrical 
waves in 3D). Similarly as before, a converging wave (cre- 
ated e.g. by the two-dimensional Maxwell's fish eye [5]) 
turns into a diverging wave at the focal point. The con- 
verging and diverging waves described by Hankcl func- 
tions are singular at the origin, but the combined wave 
satisfies the wave equation at the origin and is described 
by the Bessel function of the first kind Jo (Ax). Now, 
similarly as in the 3D case, we compare this wave with 
a plane wave of the same frequency. Fig. 2 shows the 
two waves along an x-axis passing through the focal 
point. The phase of the plane wave was chosen to co- 
incide with the phase of the cylindrical wave in the re- 
gion x J$ —A = — 27r/fc with the help of the asymptotic 
formula [9] for the Bessel function 



J ° {kx)K ^x C ° S ( kX -\) ■ (4) 

Looking at Fig. 2 and/or using Eq. (4) again, we see that 
in the region x A the cylindrical wave is delayed by 
7r/2 with respect to the plane wave, which is is precisely 
the Gouy phase shift in two dimensions. 

As the last thing we will demonstrate the effect of 
Gouy phase on sharply localized spherical optical pulses 
converging to a focal point. For convenience we surround 



the focal point at the origin by a spherical (in 3D) or cir- 
cular (in 2D) mirror of a unit radius. This makes the 
number of modes countable, which simplifies express- 
ing the pulses as superpositions of the monochromatic 
modes. We will also set c = 1 in Eq. (1), which can 
always be done by a suitable choice of units. 

Let us start with the 3D case. For simplicity we will 
consider scalar waves now to avoid relatively complicated 
boundary conditions for the vector potential [10]. Since 
we are interested in spherically symmetric pulses only, 
the spatial parts of the relevant modes depend only on 
the radius and are again described by the sine function as 
ipni?) = sm(fc„r)/r, where the numbers k n — nn, n€N 
represent the frequencies of the standing waves that turn 
to zero at the mirror. Now consider a pulse emitted at t = 
from the origin expressed as the following superposition 
of the modes 

2 ~ _ 

t()(r, t) = - 2, sm k n r sin k n t . (5) 

r n=l 

With the help of trigonometric identities and Fourier 
analysis, Eq. (5) can be rewritten as 

_^ oo 

i(>(r, t)=- ^{cos[fc„(r - t)] - cos[fc„(r + t)]} (6) 
_A(r-t)-A(r + t) 



where A denotes the periodic delta-function with period 
2,A( a ;) = E~-oo^-2n). 

We see from Eq. (7) that in the time interval t £ (0, 1) 
there is a delta peak propagating from the center to the 
spherical mirror at r = 1 that is reflected from the mirror 
at time t = 1, after which it propagates back to the 
center with the negative amplitude. The change of the 
sign is the result of the boundary condition at the mirror. 
After the converging pulse then reaches the center at 
t = 2, its sign is flipped again and it turns into a positive 
diverging pulse. However, this time the sign change is 
due to Gouy phase. The position of the peak can also be 
inferred directly from Eq. (6): it is the point where all 
the first or all the second cosine terms are in phase, i.e., 
when either r — t or r + t is an integer. 

In two dimensions the situation is more interesting. 
The rotationally symmetric modes are the Bessel func- 
tions Jo(k n r) and the frequencies k n are given by the 
boundary condition at the mirror, i.e., at r = 1, as 
Jo(k n ) — 0, so k n is the n-th zero of J . Consider a 
pulse 

oo 

II 

r) sin (k n t + — ^ (8) 

71=1 

Using Eq. (4) and the approximate values of k n derived 
from it, 

K ~ (n - - J 7T , (9) 
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we get after some manipulation 

^ oo 

«/;(r,i) « — V {cos[fc„(r - t)] + sin[fc„(r + t)}} , 

(10) 

which is valid if r is not too small. Let us now investi- 
gate propagation of a pulse described by Eq. (10). As we 
will see, the shape of the pulse is different at different 
time intervals, therefore we will discuss these intervals 
separately. 

For < t < 1 all the cosine terms in Eq. (10) are in 
phase at the point r = t. Their superposition therefore 
forms a peak similar to Dirac (5-function located at r = t 
that propagates from the center r — to the mirror at 
r = 1, see Fig. 3(a). 

To find the behavior of the pulse for 1 < t < 3, we 
make the substitution t = 2 + r in Eq. (10) using the 
approximation (9) to the frequencies k n . This gives 

^ oo 

ip(r, t) « —p V) {- sin[fc„(r - r)\ - cos[fc„(r + r)]} . 

(11) 

For t running from —1 to 0, which corresponds to 
1 < t < 2, all cosine terms in Eq. (11) are in phase 
at r = —t = 2 — t. They form a negative <5-like peak 
there that runs from the mirror at r = 1 back to the ori- 
gin, see Fig. 3(b). The minus sign of this reflected pulse 
is again the result of the boundary condition at the mir- 
ror. One might think that after the pulse passes through 
the origin at time t = 2, it will form again a diverging 
S-peak. However, this is incorrect due to the effect of 
Gouy phase. Eq. (11) reveals that for < r < 1 (corre- 
sponding to 2 < t < 3) there is no location r at which all 
cosine terms would have zero argument. However, now 
it is the sine terms that can be combined all in phase, 
namely at the point r = r, and their superposition forms 
a pulse localized around that point. But the shape of this 
pulse is now completely different from the pulse obtained 
previously by the cosine terms; instead of a 5-like peak 
we now have a double peak resembling the cotangent 
function, see Fig. 3(c), and also the localization of the 
pulse is much worse. This dramatic change of the pulse 
shape is caused by the Gouy phase of ir/2 that, for each 
monochromatic component of the pulse, converted a co- 
sine in Eq. (10) into a sine in Eq. (11). A similar pulse 
reshaping has been observed experimentally [11]. 

In conclusion, we have analyzed the Gouy phase for 
rotationally and spherically symmetric waves. We have 
shown that for monochromatic waves the Gouy phase 
has a simple meaning that can be visualized by compar- 
ing graphs for radially symmetric solutions of the wave 
equation with plane wave solutions. We have also shown 
that the Gouy phase has a dramatic effect on sharp rota- 
tionally symmetric pulses in 2D whose shape is changed 
completely after passing the focal point. 
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Fig. 3. Propagation of the pulse described by Eq. (8). The 
sum was truncated at n = 300 and smoothened by the 
Gaussian function exp[— 3.5 (n/300) 2 ] to eliminate rapid 
oscillations that would occur due to the truncation. The 
pulse is shown at times (a) t = 0.2 (<5-like peak running 
towards the mirror), (b) t = 1.2 (negative 5-like peak 
running back to the center), and (c) t = 2.5 (a cotangent- 
like peak running towards the mirror) . Arrows mark the 
direction of pulse propagation. The dramatic change of 
the peak shape after passing through the center is caused 
by the Gouy phase of 7r/2. 

ences Research Council. 
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